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Abstract 



This paper explores a new mechanism for B production in which a b quark 
combines with a hght parton from the hard-scattering process before hadroniz- 
ing into the B hadron. This recombination mechanism can be calculated 
within perturbative QCD up to a few nonperturbative constants. Though 
suppressed at large transverse momentum by a factor AqcD^Tife/p^ relative 
to b fragmentation production, it can be important at large rapidities. A 
signature for this heavy-quark recombination mechanism in pp colliders is 
the presence of rapidity asymmetries in B cross sections. Given reasonable 
assumptions about the size of nonperturbative parameters entering the calcu- 
lation, we find that the asymmetries are only significant for rapidities larger 
than those currently probed by collider experiments. 
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I. HEAVY QUARK PRODUCTION 



Perturbative QCD calculations of heavy hadron production are based on the factor- 
ization theorems of QCD H]. These theorems imply that, for sufficiently large transverse 
momentum, production of a heavy-light hadron occurs in two stages: first, a heavy quark 
is produced in a short-distance partonic process; second, the heavy quark fragments into 
the hadron. They allow cross sections to be factorized into partonic cross sections and 
fragmentation functions. For example, the formula for the B production cross section in a 
gluon-gluon collision is 

da[gg^B + X] = da[gg^bb]0Dj;^s. (1) 

(We denote hadrons with b quarks as B and hadrons with b antiquarks as B.) The symbol 
® represents convolution over the final-state momentum fraction, z, of the B hadron in 
the b quark jet. The fragmentation function Di_^^{z) is a universal function that involves 
long-distance physics and hence cannot be computed perturbatively. On the other hand, 
the partonic cross section da[gg bb] depends only on short- distance physics and can be 
calculated in perturbation theory. To obtain cross sections for hadronic collisions, partonic 
cross sections such as Eq. (|I|) are convolved with parton distribution functions for the initial- 
state partons. 

The QCD factorization theorems are exact only in the limit of large transverse momen- 
tum, p_L — >• oo. For finite p_L, one expects corrections that are suppressed by powers of 
Aqcd/p±- In this paper we will argue that there are 0{KqcY)mh/p\) corrections that are 
calculable within perturbative QCD up to a few multiplicative constants. These correc- 
tions come from processes in which a light quark, g, participating in the hard-scattering 
process combines with the b and they subsequently hadronize into a state including the B 
hadron. We will refer to this as a bq recombination process. In pp colliders, this process 
can be distinguished qualitatively from b fragmentation because it gives interesting rapidity 
asymmetries in B hadron production. Though not important in the central rapidity re- 
gion, bq recombination can actually give a larger contribution than b fragmentation at large 
rapidities. 

Before describing this mechanism in detail, it is worthwhile to briefiy summarize the 
current status of B production. Theoretical calculations of heavy quark production have 
been performed to NLO in perturbative QCD Comparison of these perturbative QCD 

calculations with collider data is reviewed in At the Tevatron ||7HTl|, the shape of 

the p±_ distribution is in agreement with experiments, while its normalization appears to be 
somewhat underpredicted by the theory. However, even at NLO the theoretical prediction 
suffers from considerable uncertainty due to dependence on the renormalization scale. Set- 
ting the renormalization scale, /zr, equal to {p\ + m^)^/^ gives a prediction that is roughly 
a factor of 2 below experiment, while for ^ir ^ (p^ + m^Y^'^ /2 the theoretical prediction 
can accomodate the data. The strong sensitivity to the renormalization scale indicates that 
higher order perturbative corrections are important. Recent studies that attempt to resum 
logarithmically enhanced higher order corrections can be found in |1l2|-|l5[|. The resumma- 



tion of higher order terms of the form a" ln"(p_|_/m;,) and a" ln"^^(p_|_/mfe) gives results with 
reduced scale dependence and the prediction of the resummed calculation is similar to that 
of the fixed NLO calculation with the lower values of /i/j. Thus it seems at least possible that 



higher order perturbative QCD calculations would account for the shape and normalization 
of the experimentally observed pj_ distribution. 

Most measurements of heavy quark production focus on the central rapidity region. Less 
is known about the cross section in the forward regions. The D0 collaboration has 



measured the production of muons coming from b decay in the rapidity range 2.4 < \y^\ < 
3.2. They found that the excess over the theoretical prediction in the forward region is 
greater than in the central region, indicating perhaps that large rapidity b production is 
not as well understood. If one chooses fin so that the central region is described correctly, 
the NLO calculation of the cross section at large rapidity is too small by a factor of 2. 
It has been suggested that modification of the heavy quark fragmentation function or 
a reorganization of the perturbative calculation [|l^] could modify the forward-to-central 
rapidity ratio. However, CDF measurements of bb rapidity correlations agree well with 
NLO QCD calculations. CDF measured the ratio of central {\y\ < 0.6) to forward (2.0 < 
\y\ < 2.6) b production subject to the constraint that the other b was produced with \y\ < 1.5 
and that both 6's had p± > 25 GeV. It is important to point out that one cannot directly 
compare the D0 and CDF experiments because a) they involve slightly different rapidity 
regions, b) D0 makes no cuts on the second b produced in the event and c) most importantly, 
the D0 experiment uses a much smaller p± cut of 5 GeV. 



II. HEAVY-QUARK RECOMBINATION MECHANISM 

Our investigation of novel mechanisms for B production is motivated in part by the 
possibility that there are important contributions to B production at large rapidity that 
are not accounted for in Eq. ([l|). This paper will consider contributions to B production 
in which a light quark involved in the short-distance process combines with the b and they 
subsequently hadronize into the final-state B meson. This bq recombination contribution 
is distinct from what is included in standard fragmentation calculations, where the parton 
process qg b + b + q is taken into account at next-to-leading order. Over most of the 
phase space, the q in the final state produces a jet that is distinct from those produced by 
the b and b. There is a small contribution from the corner of phase space where the q has 
small momentum in the rest frame of the b. However, in this region, the q and b can bind 
and large nonperturbative effects can enhance the cross section. Our mechanism takes these 
nonperturbative effects into account. 

The word "recombination" in high energy physics dates back to the parton recombination 
model invented by Das and Hwa in 1977 to describe the production of mesons at low p± in 
hadronic collisions |T^ . In this model, most of the momentum of the meson comes from the 



momentum of a valence parton in one of the colliding hadrons. The meson is formed by that 
valence parton recombining with a sea parton from that same hadron. Our recombination 
mechanism for heavy meson production is analogous in that it also allows a valence parton 
to deliver most of its momentum to the heavy meson, although in this case the valence 
parton combines with a heavy quark created in the hard scattering. Our recombination 
mechanism is also analogous to a "higher twist" mechanism for the inclusive production of 
light mesons at large p± studied by Berger, Gottschalk, and Sivers [|T^]. This mechanism 
also involves a parton from the colliding hadron that is involved in the hard scattering and 
ends up as a constituent of the large-p_L meson. The cross section is suppressed by ^qcd/ Pi. 
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relative to the fragmentation mechanism, and the nonperturbative factor involves a function 
that describes the distribution of the light-cone momentum in the meson. In our mechanism 
for heavy meson production, the cross section is suppressed only by Aqco'^fe/pi, and the 
nonperturbative factor can be reduced to a few numbers. For p± ~ rrih, the contribution 
from our mechanism is suppressed only by Aqcd/?^?)- In this region of p±, there are other 
contributions suppressed by Aqcd/w^?) in which the parton that combines with the heavy 
quark is a spectator not involved in the hard scattering []T^. String fragmentation models 



for hadronization provide an explicit model for these contributions [g^] . Another mechanism 
that gives contributions at p± ^ mh that are suppressed by Aq^j-j/m^ is intrinsic heavy flavor 
||2T| , |22[| . These other mechanisms may give contributions to the integrated cross section that 
are comparable to that of our heavy-quark recombination mechanism, but they should be 
smaller for p± ^ nib. Even for p± ~ rrih, our mechanism may be numerically important in 
certain kinematic regions, such as large rapidity, because it allows most of the momentum 
of a valence parton to be delivered to the heavy meson. 

To motivate our recombination mechanism for heavy mesons, it is helpful to make an 
analogy with quarkonium production [p^-|27[]. Quarkonium production cross sections also 



obey a factorization formula: 

da[tj -,H + X] = Y: da[tj (QQ^ + X]{0\O^\0) , (2) 

n 

where i,j represent initial-state partons and QQ represents the heavy quark- ant iquark pair. 
The QQ are produced in the short-distance process with small relative momentum in a 
state with definite angular momentum and color quantum numbers, which are collectively 
denoted by n. The QQ pair in the state n subsequently evolves into a final hadronic state 
that includes the quarkonium, H. The probability for this is given by a vacuum matrix 
element of a local operator in NRQCD, denoted by . If the QQ produced in the short- 
distance process is a color-singlet and has the same angular momentum quantum numbers as 
the quarkonium being produced, the matrix element has a physical interpretation in terms 
of the nonrelativistic wave function of the QQ pair in the bound state. 

The quarkonium factorization formula suggests that there are analogous contributions 
to B production: 

da[qg ^ 5 + X] = J] da[qg ^ (5g)" + b] p[(6g)" ^ B] . (3) 

n 

(A similar mechanism for B production is obtained by replacing b and q q). The 
notation means that the q has small momentum in the b rest frame and that the b and 
the q are in a state with definite color and angular momentum quantum numbers specified 
by n. The factor p[(bq)^ — > B] represents the probability for to evolve into a final state 
that includes the B hadron. We refer to the process in Eq. @ as bq recombination. 

Note that the color and spin quantum numbers of the (bq)"' need not coincide with those 
in the quark model for the B hadron. The (bq)"' can evolve into a state that contains the 
B hadron and an arbitrary number of partons that have small momentum in the B hadron 
rest frame. In this nonperturbative transition, soft gluons can radiate color and flip the 
spin of the light quark so the (bq)" produced in the short- distance process need not be a 
color-singlet and its angular momentum quantum numbers can differ from those of the final- 
state B hadron. In quarkonia production, these nonperturbative transitions are governed by 
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the multipole expansion and the nonperturbative matrix elements appearing in Eq. (0) scale 
with definite powers of the typical velocity, f , of the heavy quarks in the bound state. Those 
transitions for which the color and angular momentum quantum numbers of the (QQ)"' do 
not match those of H are suppressed by powers of v, and are therefore controlled by the 
heavy quark mass. We expect no such suppression for the parameters p[(bq)"' —>■ B] in 
Eq. (I). 

Furthermore, it is important to observe that the B meson need not have the same light 
quark fiavor quantum numbers as the pair that is produced in the short- distance 

process. The hadronization process can involve the creation of a light quark-antiquark pair, 
q'q', and the b can bind with the q' to form a B meson with a different light quark flavor 
from the that is produced in the hard scattering. This is an important difference 

from quarkonium production, where the total heavy quark and heavy antiquark numbers 
are separately conserved because the hadronization process is governed by long-distance 
physics. Light quark-antiquark pair production could of course have a nonperturbative 
suppression factor. For instance, there is such a suppression factor in the large Nc limit 
of QCD. Thus a (5m)" could prefer to hadronize into a B~^ rather than a B^. If such a 
suppression factor exists, it is not controlled by the heavy quark mass and should therefore 
be roughly the same for b and c hadrons. The flavor dependence of the nonperturbative 
factor p[(bq)"' B] is important for predictions of flavor-speciflc B hadron production cross 
sections. In this paper, we will for simplicity consider only cross sections that are summed 
over all B hadrons. For this inclusive cross section, the nonperturbative factor in Eq. 
can therefore be denoted simply by pn- 

It is only necessary to consider with relative orbital angular momentum L = 0. 

We will show below that the recombination constants p[{bq)" — ^ B] for L = are propor- 
tional to Aqcd/ nil,, so the recombination contributions to the cross section are suppressed by 
AqcD'^fe/pl- Producing (bq)^ in a higher orbital angular momentum state requires keeping 
additional factors of the light quark momentum in the hard-scattering matrix element. The 
light quark momentum in the bound state is O(Aqcd) and these additional factors of Aqcd 
must be compensated by powers of mi, or p±. So up to corrections suppressed by additional 
powers of Aqcri/^b or Aqcd/p±5 it is sufficient to compute production of in the spec- 
troscopic states ^Sq^\ ^S[^\ ^Sq^^ and ^S[^\ where the superscript in parantheses indicates 
either color-singlet or color-octet. Furthermore, heavy quark spin symmetry p8| , p9| relates 
the transition amplitudes of ^5*0 and ^Si states. Thus the calculation involves only two 
nonperturbative parameters, pi for color-singlet and pg for color-octet production. Heavy 
quark fiavor symmetry can also be used to relate the production parameters for b hadrons 
to those for c hadrons. 

In most regions of phase space, b fragmentation dominates over the bq recombination 
mechanism. However, the bq recombination cross section is peaked in the forward region 
where the (bq) pair emerges in approximately the same direction as the initial-state light 
quark. As shown below, it is possible for bq recombination to be comparable to b fragmenta- 
tion in the forward region. In addition to enhancing the cross section, recombination leads 
to an interesting rapidity asymmetry in B hadron production at pp colliders. Let the proton 
direction define positive rapidity. If the light quark in the initial state is from the proton, 
the (bq), and therefore the B hadron, will emerge preferentially with large positive rapidity, 
while if the initial-state light quark comes from the antiproton, the B hadron will prefer to 
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FIG. 1. Diagrams for bq recombination into a B meson. Light quarks are single lines, heavy 
quarks are double lines and the shaded blob is the B meson. 



have negative rapidity. Since the proton has more u, d than m, d, there should be an excess 
of (hu) and (M) at positive rapidity, while the u, d in the antiproton should lead to an excess 
of {hu) and {hd) at negative rapidity. This rapidity asymmetry is a striking prediction of 
this mechanism. 

Rapidity asymmetries have also been predicted on the basis of string fragmentation 
models for hadronization In these models, the asymmetry arises from a process 



called the "beam drag effect" in which the h produced by the hard scattering combines 
with a spectator quark from an initial-state hadron. Since the beam drag effect necessarily 
involves soft momentum scales, this nonperturbative contribution is part of the corrections 
to Eq. (|1|) suppressed by at least Aqqj-,/p^ at large p_L. We emphasize that our recombination 
mechanism generates the asymmetry by a perturbative mechanism in which a h combines 
with a light quark participating in the hard scattering process. Thus, it is suppressed at 
large only by Aqcd'^6/Pj_- Our process is therefore more important than the beam drag 
effect for ^ '^b- Even for px ~ m^, our process may be numerically more important at 
large rapidity, because it allows most of the momentum of a valence parton to be delivered 
to the heavy meson. 



III. HEAVY-QUARK RECOMBINATION CROSS SECTION 

This paper presents calculations of the hq recombination contribution to the inclusive 
cross section for B hadrons at the Tevatron. We will focus on the rapidity dependence of B 
hadron production. The rapidity dependence of the inclusive cross section for B hadrons can 
be obtained by replacing y —y. The leading order diagrams contributing to qg — > (bq) +b 
are shown in Fig. |I]a-e. The single lines represent light quarks and double lines are heavy 
quarks. The shaded blob represents the transition amplitude for the (bq) to hadronize into 
a state including a B hadron. 

Let us denote the momentum of the initial-state light quark as p, the initial-state gluon 
as /, the final-state light quark as Pg, the momentum of the b (which ends up in the B) as 
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Pb, and the momentum of the b quark as k. In the rest frame of the B hadron, pt has a large 
component of 0{miy) while all components oipq are O(Aqcd)- Therefore, in the perturbative 
amplitude we can take — up to corrections of ©(Aqco/'^^f))- The 0{p~^) piece of the 
amplitude is: 



M 




2p-Pg 

v{Ph)u{Pq) 




(4) 



where Ai and A2 are given by: 



= T^Yu{p) ■ u{k) (^"T^ 



- /+ ^- nib 
abcrrc ( /- :^)e/. + 2p ■ e 7^ - 2 /^ 



/- 



rUb 



(5) 



(k + Pb)' 



{k + PbY 
■u{p) ■ u{k)T'''yf, 



The trace in Eq. (Q) is over both Dirac and color indices. The polarization vector of the 
gluon is e^. Contributions to Ai come from diagrams in Figs. |l]a)-c), while A2 comes from 
Fig. Ilje). Note that Fig. |l|d) does not have any terms that are singular as Pq 0, so we 
neglect its contribution relative to the terms kept in Eq. (^. 

At this point we need to project the partonic amplitude onto a matrix element that gives 
the transition amplitude to a state with a B hadron. Let us assume for the moment that the 
(bq) evolves into the exclusive state consisting of a single B meson, which is possible only if 
the (bq) are in a color-singlet state. (Later, we will consider the contribution of (bq) evolving 
into final states containing additional soft hadrons.) Then we should replace the spinors in 
Eq. (^) with the valence light-cone wavefunction of the B meson. The formalism we use is 
similar to that used for exclusive B decays and form factors |3O-0^. We will employ the 



conventions of [32,33]. The light-cone wavefunction of the B meson is: 



[B{v)\qp,,iz)b^,iO)\0) 



5ij ifsmB 

Nr 4 



1- 



20+(t)- / } 75 



(6) 



J a(3 



In Eq. (^, z is a light-like vector, pb = rubV and t = v ■ z. The a, 13 are Dirac indices, i.j are 
color indices. For (bq) in a '^S^^ state, we should replace the B meson with a B* meson and 
the 75 with where is the polarization vector of the B* . The right hand side is the most 
general Lorentz decomposition of the non-local matrix element consistent with heavy quark 
symmetry. The functions 0±(t) are defined so that 0±(O) = 1 and fs is the B meson decay 
constant. The amplitude for producing the B meson is obtained by making the following 
substitution in the matrix element for the partonic cross section in Eq. ||33[| : 



Vi{Pb)Uj(Pq 



(7) 



rrib 



4 



dco{2M^)+ / rfr7(0_(r7)-0+(r/))7^ 



dpq 



75 



Pq=UJV 
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The functions (f){uj) are the Fourier transforms of (f){t). When making this substitution, pq 
should be replaced with uv only after taking derivatives on the amplitude. It is straightfor- 
ward to take derivatives and integrate by parts to obtain: 



4 - "^b)75 {\f+- 



where the f± are defined by 



4 JO uj 



f- 





21 -v 



f- 



21 -v 



(8) 



(9) 



id 



In the last equation, we have used 
of/- 



is divergent [32,33 



(0) = 0. Since 0-(O) ^ 0, the integral in the definition 
This divergence is cut off by the transverse momentum of the 
light quark in the B meson. If the cross section depended on /_ it would indicate that the 
recombination mechanism is sensitive to soft physics. Fortunately, this is not the case. It is 
possible to show that ^Ai = /A2 = 0. For example. 



/A2 ~ / / / = 2e ■ / /- 



0. 



(10) 



Thus all terms proportional to /_ give a vanishing contribution to the amplitude. The final 
result for the amplitude for the exclusive recombination production of a i? meson is: 



M 



B 



Tr 



4 - rnb)j5 



Ai 



+ 



Ao 



2p -pb 21- pb 



(11) 



This amplitude can be obtained from the partonic amplitude of Eq. (^) by making a few 
simple substitutions. For the Dirac spinors, we make the substitution 



Vi{Pb)Uj{Pq) 



-^mbf+ifi - rribh^ 



(12) 



and then we replace all other occurrences of Pq in Eq. (^) with pb- Since fs scales like 



^qcb/^I^"^ and (f)+{uj) scales like Aq^^, we conclude that /+ scales like (AQco/'^^fe)"'^''^- 

Note that the final expression obtained for the production matrix element is very similar 
to what would be obtained in a quarkonium production calculation. If we replaced the light 
quark with a charm quark in the above calculation, then the outgoing momentum would be 
Pq = rpB and pb = {1 — r)pB, where r = rric/ {fric + mb) and Pb is the B^. momentum. The 
projection operator for the B^. is 



v{pb)iUj{pq 



(13) 



Thus the matrix element for B production via (bq) can be obtained by taking 

the r — > limit of the matrix element for B^. production and replacing fB^/{4:r) —>■ f+mb- 
The resulting cross sections are equivalent to modelling the B meson as a nonrelativistic 
bound state with a light constituent quark. 
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In calculating the inclusive cross section, one must also include processes in which the (bq) 
hadronizes into multiparticle final states, such as Btt, Btttt, etc. After squaring the matrix 
element, the exclusive cross section is proportional to The field theoretic expression for 
fl /+ involves a projection operator | B{v)){B{v) \ . The inclusive cross section will instead 
involve the inclusive projection operator 



Y.\B{v) + X){B{v)+X\ 



(14) 



X 



where the sum is over all possible soft hadrons X in the rest frame of the B. After ma- 
nipulations analogous to the exclusive case, we would ultimately arrive at an expression 
proportional to a double moment of some density matrix instead of The bottom 

line is that the inclusive cross section for color-singlet hq recombination can be obtained from 
the exclusive cross section simply by replacing by a different nonperturbative parameter 
Pi- The parameter pi is proportional to the probability for the color-singlet to evolve 
into any state containing the B hadron. Like the parameter pi scales like AQCD/''^b- 

If the vacuum saturation approximation is applied to the sum over soft hadrons in 
Eq. ([I^) , it would give as an estimate for pi. In quarkonium production [Q, the vacuum 
saturation approximation can be applied to color-singlet matrix elements with the appropri- 
ate angular momentum quantum numbers. It is a controlled approximation with corrections 
that scale as f^. In the case of B hadrons, the vacuum saturation approximation is not a 
controlled approximation and it is likely to drastically underestimate the value of pi. There 
is no reason to expect the contributions from Svr, Stttt, etc. to be suppressed relative to B. 
The extraction of the parameter /_|_ from exclusive B decays could be used to place a lower 
bound on pi, but we emphasize that it should be regarded as a weak lower bound. 

Our final result for the color-singlet hq recombination contributions to the inclusive cross 
section for B hadron production in qg collisions is given by 



da 
~dt 



243 53 



27r^«3 ^2 



64 f/ 



S 



T 



79 



16 S 



112 U 64f/2 



T 




8U 



64 



2^2 



2^2 



48 mis 



(15) 



(16) 



8U 



where we have defined S = {p + /)^, T = t — mf = {p — PbY 
- Pbf - ml = -21 ■ pb. 



'2p ■ Pb and U 



To understand the behavior of the cross section, it is useful to 



express the variables S, T and U in terms of p± and Ay 
rapidities of the b and b, respectively: 



y — y, where y and y are the 



5" = 2{p']_ + ml){l + coshAy) 
T = -(pi + m^)(l + e-^^), 
U = -{pl + ml){l + e'^y). 
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At fixed A?/, da/dt = pi da/dt ~ AqcBmb/ {p^ + ^lY- For large pj_, the cross sections da/dt 
for recombination scale like Aqcd^6/Pj_) while the partonic cross sections that appear in 
the h fragmentation formula in Eq. (|1|) scale like Thus, recombination is suppressed at 

large p± by Aqco'^fe/pi relative to h fragmentation, consistent with the QCD factorization 
theorems. Although the recombination cross sections are peaked in the forward region, they 
are always finite as the h quark mass cuts off all infrared divergences. 

It is also possible for a color-octet to hadronize into a state with a B plus soft hadrons. 
To obtain the color-octet production cross section starting from the partonic amplitude in 
Eq. (Hp, we first make the substitution 



Vi{pb)Uj{p. 




(17) 



and replace all other occurrences of Pq with pb. After squaring the matrix element and 
summing over colors in the final state, (/+)^ is replaced with pg, a factor that is proportional 
to the probability for a (bq) in a color-octet state to hadronize into a state including the B 
hadron. The normalization in Eq. ([T7|) is chosen so that ps/pi is the ratio of the probabilities 
for a color-octet (bq) and a color-singlet (bq) to hadronize into a B hadron. Since both pi 
and ps are 0(AQCD/'"^b), we might expect this ratio to be approximately 1. Our results for 
the color-octet bq recombination contributions to the inclusive cross section for B hadron 
production in qg collisions are: 



An'^al ml 



53 




U f 9T 9T^' 

"5 v^ir^ip 



(18) 



UU 8 f/2 
79 + —— + 



T 



2^2 



|l5<;('sf')l 



243 5^ 



2T 



233 + 



U / 9T 



lOU 316 T 266 T2 8U^' 



8rn 


mi S 


f/2 ) 


+ 


18U 
T ' 


h 

f/2 



9T 8t/ 



2^2 



(19) 



+ 



[/2 



2'2 



+ 



3 S' 



9T 8U 



The normalization of the bq recombination contributions to the inclusive cross section 
for B hadron production depend on the parameters pi and pg. Since their values are not 
known, we will focus on the qualitative aspects of the rapidity distributions in pp collisions. 
Note that bq recombination can actually contribute to B hadron production in two ways: 



a) Yl M<19 ^ {bqT + b] Pn[{W ^ B] 

n 

b) 1^9 ^ (6g) " + &] p„ ® D-,^j, . 



(20) 



q,n 



In a), the {bq) hadronizes into the B hadron. In b), the {bq) hadronizes into a B, and the B 



comes from the fragmentation of the associated b. In pp collisions, both processes in Eq. (20) 
lead to asymmetric B hadron rapidity distributions. As discussed earlier, process a) yields B 
hadrons strongly peaked in the proton direction. The most important contribution to process 
b) comes when the initial-state g is a m or d from the antiproton, and the resulting B hadron 
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FIG. 2. Color-singlet (solid) and color-octet (dashed) contributions to the rapidity distribution 
of B hadrons in pp collisions at ^/s = 1.8 TeV, with a p_l > 5GeV cut. We have chosen pi = 0.1 
and ps = 0.17, so color-singlet and color-octet cross sections are roughly equal at central rapidity. 

cross section is peaked at negative rapidity. Since a) produces -B's with positive rapidity 
while b) produces 5's with negative rapidity, the asymmetry in the rapidity distribution 
tends to get washed out when the two effects are combined. The extent to which this happens 
depends on the color state of the (bq) and (bq). This is illustrated in Fig. |^, which shows the 
rapidity distribution of B hadrons from bq recombination in pp collisions at y/s = 1.8 TeV, 
with a pj_ > 5 GeV cut. We have set rrib = 4.75 GeV and used the CTEQ5L parton 
distribution functions. We used the one- loop running formula for asifiR) with parameter 
A^^) = 0.146 GeV and set the renormalization and factorization scales to (m^ +p\Y^^. The 
color-singlet and color-octet bq recombination cross sections are the solid and dashed lines, 
respectively. The unknown parameters pi and pg are 0{AQCD/mh) and therefore are expected 
to be less than 1. We have chosen pi = 0.1 and pg = 0.17. Their relative normalization is 
such that the two distributions are approximately equal at central rapidity. Color-singlet 
recombination produces an asymmetric rapidity distribution with a small peak in the forward 
direction from process a) and a large peak in the backward direction from process b). In 
color-octet production, the processes a) and b) combine to yield a surprisingly symmetric 
rapidity distribution. Thus the B rapidity asymmetries can give important information on 
the relative importance of color-octet and color-singlet production mechanisms. 
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FIG. 3. Rapidity distribution of B hadrons in pp collisions at y/s = 1.8 TeV, with a > 5 GeV 
cut. Solid line is NLO b fragmentation contribution, dotted and dashed lines are color-singlet and 
color-octet contributions. The parameters pi and pg are the same as in Fig. ^. 

IV. COMPARISON WITH FRAGMENTATION MECHANISM 

Next we compare the bq recombination cross sections with the NLO perturbative QCD 
calculations of b fragmentation production of B hadrons in pp colliders. The rapidity dis- 
tributions are compared in Fig. ^ The solid line shows the b fragmentation contribution, 
which is obtained by using the LO cross section for b production convolved with CTEQ5L 



parton distribution functions then multiplied by a K factor of 2.0 to give a normaliza- 
tion consistent with the NLO calculation In PJTO[], it is demonstrated that the shapes 
of the NLO and LO rapidity distributions are nearly identical, so multiplying the LO cross 
section by a K factor gives an adequate representation of the NLO calculation. The partonic 
cross section should be convolved with a b quark fragmentation function. We have checked 



using the Peterson fragmentation function that the shape of the rapidity distribution is 



essentially unaffected by this convolution. Also depicted in Fig. ^ are the bq recombination 
cross sections for B hadrons, including both processes in Eq. pO|). The dotted and dashed 
lines show the rapidity distributions for color-singlet and color-octet bq recombination, re- 
spectively. The parameters pi and ps are the same as in Fig. |^. We emphasize that the 
normalization of the recombination contributions in Fig. ^ is uncertain, and that only their 
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FIG. 4. Rapidity asymmetry for B hadrons in pp collisions at y/s = 1.8 TeV, with a > 5 GeV 
cut. The solid line is color-singlet, dashed line is color-octet. The parameters pi and are the 
same as in Fig. 0. 



shapes are determined. It is clear that the recombination cross section can be as large as the 
fragmentation contribution for sufficiently large rapidity, however the point at which this 
occurs depends sensitively on pi and pg- 

It is also interesting to look at the rapidity asymmetry of B production, which is defined 

by 

This is shown in Fig. ^. Color-singlet recombination, shown by the solid line, only gives a 
noticeable asymmetry for \y\ > 3, while color-octet recombination, shown by the dashed line, 
gives a negligible asymmetry for \y\ < 5. For the values of pi and ps chosen in this paper, 
the asymmetry is unobservable at the rapidities currently probed by the D0 experiment. 
The asymmetries can of course be larger or smaller depending on the values of pi and pg. 

In order to assess how important the bq recombination mechanism is in the rapidity 
regions that are accessible experimentally, one needs to know the normalization of the cross 
section, which is determined by the parameters pi and pg. Unfortunately, these parameters 
are unknown. The vacuum saturation approximation gives the weak lower bound pi > fl. 



13 



If one makes the mild assumption that the hght-cone wavefunction (f)+{uj) is positive definite, 
it is possible to place a lower bound : /+ > 0.1 |31].[| The light-cone wavefunction models 



employed in [^] give estimates for /_|_ in the range 0.1 — 0.2. We conclude that pi > 0.01. 

The actual determination of pi and pg would require phenomenological analysis of B 
hadron production in other high energy physics processes. We can extract estimates of 
these parameters from a constituent quark model for heavy-meson production developed 
by Berezhnoy, Kiselev and Likhoded PB|-P^. In this model, the B hadron cross section 
is obtained from a factorization formula analogous to Eq. (|^) for be hadrons by replacing 
rric with a constituent quark mass, rriq, for the light quark and computing the parton cross 
section to leading order in perturbative QCD. In the simplest form of this model, only 
the color-singlet S-wave matrix element, (Of), is kept and it is estimated from quark model 
wave-functions [Q. This model was previously used to calculate heavy quark fragmentation 
functions . However, the model can easily be extended by relaxing the vacuum saturation 



approximation for (Of) and including a color-octet S-wave matrix element (Og) |37-39 



The parton cross sections in this model can include fragmentation, recombination, and other 
contributions, depending on the high energy process being considered. All the contributions 
are calculated in terms of the two matrix elements (0f ) and (O^). Our recombination cross 
sections can be obtained from calculations using this model by isolating the recombination 
contributions and then keeping the leading terms in the limit r = mq/{mq + rub) — > 0. We 
can therefore use the matrix elements of |37|^9i to estimate our parameters pi and ps- 



The simplest way to determine the values of pi and ps that correspond to the matrix 
elements (Of) and (Og) of |l37|^9| is to compare their calculations of D* fragmentation 



functions with what we would get if we applied our formalism. For the parameters 
Pn[cq — > D*], we find pf^g = 3(Oi^g)/(4mcmg), where nig = 0.3 GeV is a constituent 
quark mass. /^From studies of photoproduction and electroproduction PB| , |5D| , they extract 
(Ci) + (Cg)/8 = 0.25 GeV^ and (Og) = 0.33-0.49 GeV^. The corresponding values of the 
parameters in our formalism are pf -|- pf /8 = lA,pg* = 1.8 — 2.8. We can estimate the 
corresponding parameters pfg for B* production by multiplying by a factor of mc/mh. This 
is because p^g scale with the heavy quark mass as l/rrib- Note that we are assuming that 
the moment of the light-cone wavefunction appearing in the definition of /+ does not exhibit 
a strong dependence on the heavy quark mass between the scale rric and mf,. The resulting 
estimates are pi ~ 0.4 and pg ~ 0.8, which are larger than the values we used in the figures. 

So far we have focused exclusively on the contributions in which the parton that combines 
with the b is the light quark q. There are also processes in which a q participates in 
recombination: 

da[qg ^B + X] = J2 d(j[qg ^ + b] pn[bq ^ B] . (22) 

n 

The LO diagrams for the process in Eq. ( p2|) are identical to the diagrams shown in Fig. |l|, 
except the light quark lines are replaced with light antiquark lines. Since the diagrams for 
bq recombination are almost identical to those for bq recombination, we expect that the 
cross section will be peaked in the direction of the incoming q. The B hadron will emerge 



^We thank Dan Pirjol for discussions on this point. 
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preferentially in the p direction rather than the p direction, so one would expect that the 
asymmetry would be the reverse of that seen in Fig. ^. In the case of color-octet production, 
we saw that the asymmetry was almost completely washed out by the process b) in Eq. (pOD. 
Thus, it is difficult to tell how big the asymmetries from bq recombination will be without an 
explicit calculation. The bq recombination process will introduce additional nonperturbative 
parameters. To lowest order in AqcB/^b and Aqcd/p±, the b and q are in an L = state. 
Their color state can be either the 6 or the 3 of SU(3). Heavy quark spin symmetry relates 
the spin states so there will be two parameters, ps and p^, that are analogous to the pi and 
Ps parameters for bq recombination. Intuitively, one might expect (bq) to prefer to hadronize 
to B mesons and the diquark (bq) to prefer to hadronize to the A;, by picking up another 
light antiquark. This intuition leads one to expect that bq recombination is more important 
for heavy meson production and that bq recombination is more important for heavy baryon 
production. This would probably lead to rapidity asymmetries that have opposite signs for 
Afo baryons and B mesons. 

One could also consider processes in which a gluon from the hard scattering combines with 
the b. These do not produce any interesting asymmetries because the gluon distributions in 
the p and p are identical. The (bg) can be in a 3, 6 or 15 of SU{3), so three nonperturbative 
parameters are required to describe this recombination mechanism. At large rapidity, the 
parton cross sections for (bg) may be larger than those for (bq). For reasonable values of pi 
and ps, the bq recombination cross section is too small to explain the excess b production 



at large rapidity observed by the D0 collaboration [|rT|. However, it is possible that the 
discrepancy could be explained by bg recombination. 

The excess b production at large rapidity observed by the D0 collaboration cannot be 
explained by the conventional fragmentation mechanism. We have proposed a perturba- 
tive recombination mechanism that is suppressed by Aqco^b/pi at central rapidity, but 
may give an important contribution to the cross section at large rapidity. In pp collisions, 
our mechanism generates asymmetries in the rapidity distributions of B hadrons. The ob- 
servation of rapidity asymmetries at the Tevatron that are too large to be explained by 
string fragmentation models of hadronization would provide qualitative evidence for our 
mechanism. Our mechanism involves several nonperturbative parameters, so quantitative 
evidence would require comparing its predictions with data on B production in other high 
energy physics experiments. We have calculated the cross sections for bq recombination in 
qg collisions and studied their effects on B production at the Tevatron. Future work should 
involve a complete calculation of the other recombination contributions discussed in this 
paper, and an exploration of the implications of heavy-quark recombination mechanisms for 
b and c hadron production in ep collisions and fixed target experiments. 

We thank S. Fleming, D. Pirjol and D. Soper for useful discussions. This work was 
supported by the National Science Foundation under Grant No. PHY-9800964 and by the 
Department of Energy under grant DE-FG02-91-ER4069. 
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